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Abstract 

We cafcufate the opticaf diffraction radiation generated by a bunch of high energy par- 
ticies as they pass through a round hofe within an annuiar metalfic ring. We derive ex- 
pressions for the differential angular spectrum in the far-field and the intensities of the 
horizontal and vertical poiarizations. The sensitivity of the spectrum to changes in beam 
size and position is shown. The totai photon yieid from the bunch is caicufated and used to 
set fimits on the detectabie wavefengths. 
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1 Introduction 



The use of optical diffraction radiation (ODR) as a diagnostic tool has increased in recent years. 
The potential of this technique has been demonstrated in several experiments at KEK [T], APS 
[|2J, FLASH [3J and possibly other facilities. These experiments were performed in extraction 
beam lines of lepton machines. However this technique can also be applied to high energy 
hadron beams. In this report we consider the ODR produced by such beams with the target as a 
round hole and apply the results to the Tevatron. 

This radiation is produced when a beam passes in the vicinity of a conducting target. The 
electro-magnetic fields due to the beam induce currents on the target and as the beam propa- 
gates, the currents change in time producing radiation both in the direction of beam propagation 
and along the direction of specular reflection from the target. This latter radiation, also termed 
backward diffraction radiation (BDR), is more useful for diagnostics since it can be directed 
out at the same longitudinal location as the target. This radiation is different from optical tran- 
sition radiation (OTR) in which the beam passes through a metal target. Transition radiation 
is not suitable for continuous monitoring of a beam in a collider due to the beam energy loss 
and emittance growth and the fact that the target may be damaged. However the techniques for 
analyzing ODR are similar in many respects to those for OTR. 

Measurements of the radiation intensity either in the near field or far-field have been used 
to determine beam positions and sizes. For example, the beam size and beam position of a 1.28 
GeV electron beam were measured in an extraction beam line at KEK [IJ using the far- field 
angular distribution of the radiation. The near-field image was used to monitor the relative beam 
size of a 7 GeV electron beam in the extraction line at APS |l2l. In principle, measurements of 
the beam divergence are also possible using the interference of ODR between two targets, as 
has been done with OTR. 

This paper is motivated by the desire to use this technique in colliders, especially for the 
LHC and possibly for future colliders envisaged such as the muon collider. A brief report on 
these prospects was presented earlier If the technique yields beam measurements with suf- 
ficient accuracy and reliability then the non-invasive nature would allow continuous monitoring 
during the length of a luminosity run. This would be valuable if the beam can be imaged close 
to the interaction points. 

Synchrotron radiation is already used as a non-invasive diagnostic tool in the Tevatron and 
will also be used in the LHC. The principal advantage of ODR is that it can be generated in a 
straight section and therefore used for imaging in an experimental insertion. The disadvantage 
is that the ODR flux is less copious than synchrotron radiation (OSR) and imaging will take 
longer than with OSR. 

In Section 2 we briefly discuss the parameters of different hadron colliders. In Section 3 
we derive the basic results for the angular differential spectrum of ODR from a round hole 
due to a bunch. We apply these results in Section 4 to find the sensitivity of the spectrum to 
beam size and offset changes. In Section 5, we calculate the expected photon yield from a 
bunch per turn as a function of frequency and we use this to find the frequency range where a 
sufficiently strong ODR signal can be obtained. In Section 6 we do a brief comparison of the 
ODR spectrum with the OTR spectrum. We briefly list in Section 7 the experimental issues 
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Tevatron 






-L/iici^y Lvjc V J 


yo\j 


250 


7000 


Bunch intensity 


2.7xl0ii 


2xl0ii 


1.1x10" 


Beam size [jUm] 


400 


1012 


807 


Beam div/opening angle 


2.9x10-3 


1.2x10"^ 


5.7x10-3 


Number of bunches 


36 


55 ^ 120 


2808 


Revolution frequency [kHz] 


47.6 


78.2 


11.2 



Table 1 : Table of parameters for hadron colliders 



associated with measuring ODR when two beams are present. We end with our conclusions in 
Section 8. We will use CGS units throughout. 

2 Hadron colliders 

Optical transition radiation (OTR) has been used in the Tevatron at injection energy to image 
the beam [6J. At collision this technique is not feasible both because of the impact on beam 
quality via multiple scattering in the target and the damage to the target itself. However ODR 
is non-intercepting and has the potential to be a useful diagnostic tool at collision energy. This 
technique also has potential in the LHC where we envision that placing ODR targets on both 
sides of the interaction point (IP) and before the first interaction region quadrupole would allow 
a non-invasive measurement of the beam size at the IR At RHIC the energy is lower so one 
would have to use longer wavelength ODR for a substantial radiation flux. 

Table \T\ shows some of the key parameters for these hadron colliders. The beam size in 
the Tevatron was calculated at CO while for RHIC and the LHC, locations in front of the first 
interaction region quadrupole were chosen. The ratio of the beam divergence to the opening 
angle of the radiation (~ I/7) is very small in all the colliders, hence the distortion of the 
spectrum due to the beam divergence should be negligible. 

3 ODR from a round hole 

The fields induced by a beam as it passes through a hole depends on the beam energy, the beam 
size, the beam position relative to the center of the hole and the shape of the hole. In this paper 
we consider a round hole as a target. First we analyse the fields from a single particle and 
generalize results obtained many years ago by Ter-Mikaelian [5J. Next we consider the ODR 
fields and spectrum generated by a Gaussian bunch of particles. 

3.1 Single particle fields 

Consider the case where a single particle moving at constant velocity v goes through a round 
annulus made of conducting material with inner and outer radii of a,„ and aout respectively. 
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The fields of the particle induce fields on the surface of the annulus. We introduce the Fourier 
transform of the fields as 



Ex = J Ea,,xe-''"dco Ey= j Eco,ye-'""d(0 (1) 

The Fourier transformed transverse fields of a particle moving at constant velocity along the 
z axis are given by L5J 

Ttv p 

£a,,, = — -e'"^/"i^i(ap) (2) 

TTV p 

The origin of coordinates is at the center of the hole and in particular the hole is in the z = 
plane, q is the particle charge, and 

a = ^ = l, P = {x'+fy/' 
vy 7 

^1 is a modified Bessel function of order one. 



3.1.1 Particle at the center of target 

We will first consider the simpler case of the particle moving the center of the target. We will 
derive equations for the fields and angular spectral distribution which will serve as a useful 
check of the more general case when the particle is offset from the center of the hole. 

Consider the field at any arbitrary point on the surface of the hole. Given the axial symmetry 
of the target, we use polar coordinates. The coordinates of a point {x,y) on the target are 



x = pcos^, y = psin^ 



Then 



E(o,x 
Em,y 



qa 

Ttv 



Kliap] 



COS0 

sin^ 



(3) 



We will calculate the fields at an arbitrary location using scalar diffraction theory. Within 
this approximation of assuming scalar diffraction theory to be valid, the fields from the entire 
target at an arbitrary observation point P can be found by integrating over the annulus 



Eco,x 

Eco,y 



ik qa , /"^^ , ^ e*^' ^ 

— — / pdp / d^-—Ki{ap 

2k Ttv Jain JO R 



COS^ 

sin0 



(4) 



where R' is the distance from the point on the target to the observation point P. If {x,y,z) are the 
coordinates of P, then 



i?' = - p cos + - p sin(|))2 +^2] 1/2 



(5) 
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Far field spectrum 

The observation point P is assumed to be sufficiently far from the target so that all points on 
the target have nearly the same phase from P. In this case the linear dimensions of the target are 

small compared to the distance from the target. This is the regime of Fraunhoffer diffraction. 

If R is the distance from the center of the hole to the point P, i.e. R = [x^ +3;^ then 
we assume here that Uout -C R. Thus in the phase term e^'^^ we expand 

R' ~R[l - -^(xcos0+3;sin0)]^/^ ~i?-psin0pcos(0p-0) 



where we define 



Then 



E(o,x 
E(o,y 



Op 

X = Pp COS 0p, y = Pp sin 0p, sinQp=—, k = ksmQp 

R 



2k %v R 



ik qoce'^^ f^""' f^^ 

— / pdp / J(^i^i(ap)exp[-j^pcos(<^p-<^)] 

R Jain 



cos ^ 
sin0 



(6) 



(V) 



To do the integrals, we use the integral representation of the integer Bessel functions 



i f 

Jn(z) = — exp[— /zcos0]cosn0J0 
2n Jo 



Then 



27t 



COS0 

sin0 



COS0P 

sin0p 



(8) 



(9) 



exp[— jlpcos((^P — 0)] 

The integral over the radius yields 

r^out _ 1 _ _ 

/ pdpK\{ap)Jx{kp) = — ^ {aout[kJ2{kaout)Kx{aaout) - aJ\{kaout)K2{aaout)\ 
Jain k^ + a-^ 

-Gin [kJ2 (kain)Ki {aatn) - aJ\ (kain)K2 (aa,„)] } (10) 

Using the recurrence relations 

xJ„+i{x) = 2nJn{x) -xJn-i{x), xKn+i{x) = 2nKn{x) -xKn-i{x) 

we can write 

fOout _ 1 

/ pdpKx{ap)Jx{kp) = ^ [T{aout;k) - T{ain;k)\ 

where 



T{a;k) = -a[kJo{ka)Ki{aa) + aJi{ka)Ko{aa)] 



(11) 
(12) 



Thus the fields are 

E(o,x 
E(a,y ^ 



kg a e 



ikR 



1 



nv R k^ + a^ 



[T{aout;k)-T{ain;k)] 



COS0P 

sin0p 



(13) 



The Poynting vector is 
—E 



^ExB* = ^[-pE,E*x-pE,E;y + P{\E,\^ + \Ey\^)z] (14) 



where we have used 

B, = -^Ey, By = l3E„ B, = 

The total energy deposited by the fields onto an element of area dA is the time integral of 
the projecteg Poynting vector 

= / dtSn (15) 

dA ,/-oo 

where h is the unit normal to the element. For an element orthogonal to the direction of propa- 
gation or direction of specular reflection 

dW f°° Be f°° f 

— = J^JtS,= !^j_Jt{\E^\^+\Ey\^) = 27tJd(o[\EcoM^ + \Eco,y\'^^ (16) 
Hence the differential angular spectrum is 

"^'^ hcR^[\Eco,,\^ + \Eco,y\^] (17) 



dOJco T 

where d^ is the solid angle subtended by the element at a distance R from the source. 
Thus in the far field, the differential spectrum is 

Define a critical frequency ©c, and dimensionless parameters M,?,g as 

«c = ^, « = ? = rsin0p, (19) 

Then other parameters can be written in terms of these dimensionless parameters as = yuj (Pain), 
a = u/ (Pain) etc and 

[T{aout;k) - T{ain;k)f = {^f i^g[tJo{^gut)Ki{jgu) +Ji{jgut)Ko{^gu)] 

-[tJo{^ut)Ki{ju) +Ji{^ut)Ko{ju) 
The angular spectral distribution thus is 

d^W ^ ^{g)^-^{g[tJodgut)KAgu)+Ji{L^^^ 



-[tJQ{jUt)Kx{ju)^h{^ut)KQ{ju)]^ (20) 
We comment on some features of this expression 
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• This spectrum depends on the magnitude of the inner radius a,„ only through the critical 
frequency (Oc. 

• The main dependence of the spectrum on the size of the target is through the dimension- 
less ratio g = aout/'^in- This is important since it suggests that the target hole may be 
enlarged to allow more space for the beam while at the same time increasing the outer 
radius without changing the spectrum. The parameter that will change in this case is the 
critical frequency COc and consequently the dimensional frequency CO. 

Th differential spectrum may be found by integrating over the solid angle 

t 



dW _ 
dco 

Define the function 



d^W An fy d^W 

-d^sinOpdOp 



dQ.dco 



72 Jo dQdco ^\-t^lyl 



dt 



(21) 



dt 



212 



v^l-f2/72 [1+r 
-\tJo{jUt)Kx{^u) +Ji{^ut)Ko{^u) 



8[tM^8i^t)Kii^gu)+Ji{^gut)KQ{^gu)] 



(22) 



The number of photons AA'^ emitted by a single charged particle into a bandwidth Aco is 



AN = — Acq = (- —)Aco 

dco hco dco 



which can be written as 




(23) 



(24) 



where Of/ = q^/{hc) ~ 1/137 is the fine structure constant. This depends on the relative fre- 
quency u = CO/cOc and the relative bandwidth Aco/co. 



Example: Tevatron parameters 

We evaluate the spectrum and the number of photons per particle for the Tevatron. Energy 
= 980 GeV, number of particles per bunch Np = 2.7 x 10^^. Figure \T\ shows the differential 
angular spectrum as a function of the parameters g, t with the frequency fixed at CO = COc- Figure 
[2] shows the spectrum as a function of the parameters u,t at a fixed ratio g = 1.5. See the figure 
captions for comments. 

Figure [3] shows the number of photons per particle calculated using the expression Eq (l24l) 
and assuming O) — COc and a 1% bandwidth or Aco/co = 0.01. The curve again shows that there 
is little gain in intensity when the target size increases beyond ^ > 2.5. After some transverse 
distance from the particle, its field has dropped to sufficiently low values that no radiators in 
the target can be excited and therefore there is no further increase in the ODR radiation with 
increasing material in the target. The number of photons emitted by a single particle in one pass 
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Figure 1 : The differential angular spectrum as a function of the ratio of the outer and inner radii 
of the target g = aout/cim and the angular variable t = ysin Op at constant u = I. Note that the 
spectrum saturates as a function of g for ^ > 3. 




Figure 2: The differential angular spectrum as a function of the ratio u = CO/cOc and the angular 
variable t = ysinOp at constant g = 1.5. The spectrum peaks close to m = 1, i.e. close to the 
critical frequency cOc- 
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Figure 3: The number of photons per particle as a function of the parameter g emitted at the 
critical frequency cOc. We observe that for g > 2, the number of photons increases slowly with 
8- 



through the center of the target can be found from this curve which is calculated for the critical 
frequency (Oc- For example at 

^=1.1 ^ AiV= 1.44x 10"^ 

A very simple estimate for the number of photons emitted by a bunch in a single pass at the 
frequency cOc is therefore 

AA^(bunch) =NpX AN = 3.2 x 10^ 

This in fact is an underestimate since it assumes that all particles are at the center of the 
target and therefore furthest from the material of the target. A more precise estimate using the 
density distribution of the bunch will be obtained in the following section. 



Near field spectrum 

Here we calculate the field distribution at a distance close enough to the target that the phase 
differences between different points on the target to the observation point is significant. This is 
the region of Fresnel diffraction. 

Here in the expansion for R' we keep the next order term in p/R. Thus 

R' = R[l - cos(<|) - 0p) + 1/2 ^ ^ _ p 0^cos((|) -^p) + \y 

The integration for the fields contains the extra phase factor exp[ikp^/ {2R)] when compared to 
the fields calculated in the far field approximation. 
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Define the dimensionless variables 



p R _ a 

P = —, r= — , a = — , 

^in ^in ^in 



aput _ 7if_ 

1 'I o o ' 

Gin 215 r 



kp^ 2 
2^ = ^^ 



(25) 



Define the complex function 



If" - kp^ 
2" pJi{kp)Ki{aa)tx^[i—]dp 



1 1 2 

pJi{-l^utp)Ki{-up)exp[irip ]dp 
op P 



(26) 



Then following similar steps as in the previous section, it follows that the Fourier transforms 
of the transverse electric fields are 



E(o,x 

Eco,y 



kqa e 



ikR 



nv R 



The angular spectral distribution is 
(fW 1 



[S{g-Xr)-S{\-Xr)] 



iyu^f\ [S{g-k)-S{\-k 



(27) 



The frequency spectrum is found by integrating over the solid angle 

dW An l-y d^W t 



(28) 



da 72 Jo dO-do) ^\ -t'^/y^ 



dt 



(29) 



Define the function 



D{g,u,r) 



dt 



Vl-f2/y2 



[S{g-Xr)-S{l-Xr 



(30) 



The number of photons emitted by a single charged particle into a bandwidth Act) is there- 



fore 




(31) 



This depends on the inner radius a,„ through the scaled variables g = aout/am,r = R/aimU = 
co/cOc- 

3.1.2 Particle offset from the center of target 

Our final aim is to find the spectral distribution from a bunch of particles. Towards that end 
we first need to know the field distribution from a particle offset from the center of the target. 
The target is in the z = plane and the particle moves with uniform velocity v along the z axis. 
Figure |4] shows a sketch of the target with center at O, the particle is at B and the field on the 
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0Q= Pp 



Figure 4: The target centered at O is an annulus of inner radius a« and outer radius Uout- B is 
the location of the particle offset by a distance b from the center of the target, A is the arbitrary 
location on the target for the field calculation. In the figure on the right, P is the point of 
observation, Q is the projection of P on to x — y plane. 
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target is calculated at point A. Again we use polar coordinates: particle at B has coordinates 
{b,x) while point A has coordinates (p,^) with respect to the center O. The distance of A 
relative to the particle at B is 



r± = [p^ + -2pbcos{x - (t>)]^/^ 
while the separation along the {x,y) axes are individually 

x = xa —xb = pcos0 —bcosx, y = yA —JB = psin^ —b^inx 
Hence the Fourier transforms of the transverse fields at A are 



qaKi{ap±) 



pcos^ —bcosx 
p sin^ —bsinx 



(32) 



(33) 



(34) 



nv p± 

As before R' is the distance of the observation point P from the location of the field, thus 

R' = [{xp-pcos^f + {yp-psm^f+zj>]^/^ 

where {xp,yp,zp) = (p/>cos0/>,p/>sin^p,z/>) are the coordinates of the point P. Figure [5] shows 
the relevant angles Op.^p. 

Integrating over the annulus, the fields from the entire target at the point P are (using scalar 
diffraction theory) 



Eo),x 


ik qa f"'>"' f 




In nv Jain Jo 



e*^'i^i(api) 



pdpd^ 



R' 



pcos^ — Z^cos;^ 
ps'm^ — bsinx 



(35) 



Far field spectrum 



We assume that the point P is sufficiently far from the target that the far field approximation 
is valid. The phase term is expanded as 



JkR' ikR 

= -ikpcos((j>-<j>p) 

R' R 



where k = k sin Op, Op is the angle made by OP with the z axis or sin Op = pp/R. 

The integrations are simplified if we write the integrands as derivatives with respect to vari- 
ables that are not integrated. We note first that 



d_ 
1 



Koiapj 



aKx{ap^ 



pcos{^-X)-b 
bp sin{^-x) 



Next we transform the fields to a frame rotated by angle X with respect to the {x,y) axes, i.e. 

(36) 



'El ' 




cos;i: 


sin;i: 




E(0,X 






-sin;i: 


cos;i: 




Eco,y 
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Then using the derivative expressions above, we find 



El 
El 



ik q e 



ikR 



Innv R 



d 

LA. 

bdx 



a out r^n 



a:„ Jo 



' dpd^ pKQ{apA_)e-'^P'°<^-^''^ (37) 



The double integral can be factorized into the product of single integrals by using the expansions 

oo 

KQ{a[p^ + b^~lbpcos{<^-X)Y'^) = £ In{ab)Kn{ap)e^<^-^^ 

n=—oo 

oo 

cxp[-ikpcos{(^-^p)] = £ {-i)V„ikp)e'<^-^''^ (38) 
The integration over is trivial and the double integral reduces to 



H 



Clin Jo 



dpd^ pKo{ap±)e 



-ikpcos{(j)—(l)p) 



2n £ {-iyin{ab)e-'"^^-1"'^ 
pJn{kp)Kn{ap)dp 

The integration over p can be done symbolically using Mathematica which yields 
pJn{kp)Kn{ap)dp = - [kJn+]_(ka)Kn{aa) - aJn{ka)Kn+i{aa)] 



k^ + a^^ 
a 



[kJn-i{ka)Kn{aa) + aJn{ka)Kn-{{(JCa)] 



where recurrence relations for the Bessel functions were used. Then defining a function r„(a; k) 
similar to the one defined in Equation (fT2l) in Section l3.1.1[ 



Tn{a\k) = —a[kJn-i{ka)Kn{(Xa) + (xJn{ka)Kn-i{(Xa)] 
After some further simplifications we can write 

n=0 

where 



72 2 CJniccb) [Tn{aout;k) - Tn{ain,k)] cosn{x - ^p) 



(39) 



(40) 



C„ = 1 for R = 0; C„ = 2(-/)" for n > 1 



and the rotated fields are 



El 
El 



ikaq 



JkR 



Cn[Tfi(^ClQiit',k^ Tn{ciifi\k^] 



I'^{ab)cos[n{x-^p)] 

{ab 
ab 



L -„^sin[n(z-<^)/> 



The fields in the lab frame (fifo xj^^coj) are obtained by applying the inverse rotation 



(41) 



E(o,x 






-sinx 




'El 


E(o,y 




s'mx 


cos;t: 




El 



(42) 
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Limit of zero offset 

Before proceeding further, we first check that the expressions derived here reduce to the 
expressions derived in the previous sub-section in the limit that the offset is zero. 
First we note that 



lim 



Eco,x 




'El ' 


Eco,y 




.^2 . 



Now using the fact that 



In{OCb) 

lim 

b^o b 



a 



Then only the n = \ term in the sum contributes and we have 



lim 



kqa 



JkR 



■[Ti{aout;k)-Ti{ain;k)] 



COS0P 



(43) 



These expressions agree with the expressions in Equation (1131) derived earlier. 

We return now to the case with offset. The differential angular spectrum is proportional to 
the square of the absolute norm which is preserved under rotations, 



co,y\ 



\Ei\^+\E2\' 



Hence the differential angular spectrum is given by 



The spectrum from a single particle therefore is 



1 



kqa 



where we have defined 



U,n„{C{ab)l'„{ab) cosm{x - ^p) cosn{x - (^p) 
I,n{ccb)In{ab) 



+mn- 



a^b^ 



■ smm{x - <^p) smn{x - (jfp)} 



It is helpful to use the scaled variables g.,u,t introduced in Equation (fT9l ). Then 



k = ^ 
Pain 

Furthermore 



k = ksinOp 



tu 



a 



kqa 



n:v{k^ + a^] 



.qr.2_}__ 



\Tffi{ciout\k^ Tm{(^in\k)] 



/3 



^gut 



-gu. 



o U[^-/m-i(V)^-(ir)+-^-(V)^--i(ir 



/3 ^"-^ /3 



/3 



/3 



(44) 
(45) 

(46) 
(47) 
(48) 

(49) 



(50) 
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Hence Umn{g,u,t) is only a function of the scaled variables g,u,t but does not depend on the 
inner radius a^. Similarly we define the scaled offset 

bs = — , ^ab = bs^ 

Clin P 

This shows that Vmn{g, u,t,b^^ is also independent of the inner aperture 

Since V„j„ is real and symmetric under the interchange of indices, we can write 

m=On=0 

Using the definition of Cm we obtain 

c +r* 



m=On=0 

m odd 
4(-l)'"/2 meven 
m — n odd 



m — n even 



Hence 



dQ.d(0 



I partic 



le{g,U,t,b,X) 



+4 I I (-1)''"-""'V„„] 

ra= 1 ;3= 1 ; | m— n | =even 



(51) 



In the limit that the offset goes to zero, this reduces to 



(52) 



which agrees with the expression in Equation (|20l) . 

The dependence of the single particle spectrum on the parameters {g,uj) is similar to that 
seen in Figures [T] and [21 Figure [6] shows the dependence of the single particle spectrum on the 
angle of observation for two different offsets b from the center of the hole. With an increased 
offset the particle is closer to the material of the target resulting in a larger radiation flux . 

The intensities of the different polarizations can be found from the components of the elec- 
tric field. From Equation l42l it follows that 

JkR 



/7o ON „ y^Cn[-^n(^OMf'^) T„{cii,-i',k)] 
f IniOCb) 

X <^ l'^{ab)cos[n{x-^p)]cosx+n-^^—-^sm[n{x-^p)]smx 



ab 



(53) 



ikaq 



JkR 



/7T ->\ r, z_^Cfi[Tn{ciout'ik) Tfi(^aifi',k)] 

C Iniocb) 
X j I'niab) cos[n{x - (I>p)] sin;t - ^ "^^ sm[n{x - ^p)] cosx 



(54) 
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-4 -3 -2 



Figure 6: The differential angular spectrum from a single particle vs the angle of observation 
{t = ysin Op) with two different offsets b = and b = 0.5a,„. Other parameter values are fixed 
at^ = 1.3,M = l,X = ^i^p = Tt/A. With non-zero offset, both the minimum and the maximum 
value of the flux increases. Note that the position of th maximum also changes with the offset. 



Hence it follows that the intensities of the two polarizations are given by 



dQdco 



,..2,4,... 
m=ln=l; \m--n\=even 



(55) 



where 



yx 

' mil 



mn 



U„,n[l'm{f^b) co^[m{z - ^p)] co^X + m ""^^!^^ sin[m(^ - ^p)] sin^] 



ab 

X [I'„{ab) cos[n(x - ^p)] cosx+n " , sin[n(;t: - 0p)] sin;i:] (56) 

oct> 

U,nn[C{^b) cos,[m{x - ^p)] sin;t - in-^^^^^^ sm[m{x - (pp)] cosx] 
X [l'„{ab) cos[n{x - ^p)] smx-n^^sm[n{x - </>p)] cos;t; (57) 

We will use the polarized intensities to examine their sensitivity to beam parameters in the 
next section. 



3.2 Spectrum from a bunch 

So far we've dealt with the spectrum from a single particle traveling through the hole. Now 
we'll consider the spectrum from a typical bunch. We assume here a Gaussian distribution of N 



16 



particles with transverse rms sizes CT^, Oy We also assume that the bunch center is offset from 
the center of the target with offsets (xQ.yo). Then 

pi^^y) = o g^p[ 



N 

= 2j^(j (j P^i^^^X:XO,s:yO,s:C>x,S:(yy,s) (59) 

The last equality defines the scaled density ps and we have scaled the other variables by a,„, the 
inner radius of the hole, 

bs = — , .^0,^ = — , yo,s = — , <yx,s = — , Oy^s = — (60) 

^in ^in ^in ^in ^in 

The differential angular spectrum averaged over the bunch distribution is 

d^W . f , f , . (fiW . 1 f. „ f . „ . d^W 



dQ.d(0 



bunch 



f f d W of f d W 



N f f d^W , 

/ bsdbs / dX Ps .^j \particle (61) 

sOy^s J J dildco 



It is important to note that the spectrum depends only on the scaled variables introduced in 
Equation (l60l) but not on the absolute values of (cr^, Oy), [xQ^yo). Hence this is a universal 
expression; the only dependence on machine specific parameters is on the beam energy and the 
bunch intensity. 

The two-dimensional integrals over {bs,x) can be factored as the product of single integrals 
over bs and x individually. However the integration over x introduces a triple summation and 
the integrals over b cannot be performed analytically. Instead we will evaluate the integrations 
over the bunch numerically. 

The bunch spectrum has the same dependence on the target size ratio g, the scaled frequency 
u and the scaled angle of observation t = ysinOp as the single particle spectrum. Figure |7] and 
[8] show the dependence of the spectrum on (gj) and (uj) respectively. As before, the intensity 
initially increases with g but flattens for g > 2.5. In the sequel I will set ^ = 1.3 in order to 
limit the size of the target. Larger sizes than this if feasible result in significantly larger ODR 
intensities mostly at low frequencies but do not change the intensity much at high frequencies. 
Figure [8] shows that the angular spectrum as a function of frequency peaks in the vicinity of 
M = 1 or CO = (Uf when viewed at the angle of maximum intensity corresponding to t = 1.6. At 
other angles, the first peak moves to other values of u and the peaks are of comparable height. 



4 Sensitivity to beam parameters 

The beam parameters that we wish to measure with ODR are the beam sizes and the beam 
positions. The angular spectral distribution is sensitive to these parameters and we examine 
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Figure 7: The differential angular spectrum from a bunch as a function of ^ = aout/ciin and the 
angular variable t = ysinOp at constant w = 1. As a function of g, the spectrum is relatively flat 
for g > 2.5. 




u 



Figure 8: Top: The differential angular spectrum from a bunch as a function of the scaled 
frequency u and the angular variable t = ysinGp at constant g = 1.3. Bottom: The spectrum 
as a function of m at f = 1.6, and other parameters at the same values. The spectrum peaks at 
around u= 1.1. 
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Figure 9: The differential angular spectrum of the horizontal polarization d^S^ /dQ.dO) for two 
different beam sizes. The red curve is obtained for the scaled beam size Ox,s = 0.1 while the 
blue curve is obtained for Ox,s = 0.12. Both the minimum and maximum values change with 
the beam size, the minimum increases more rapidly on a relative scale. 

here the dependence on these parameters. Since the detector will have a finite bandwidth in 
frequency acceptance, we include this in our analysis. We define a finite bandwidth spectrum 
by integrating over the frequency as 



Here we assume a 1% bandwidth for Au. 

Consider the sensitivity of the separate polarizations to changes in the beam size. Figure |9] 
shows the horizontal polarization intensity d^S^/dQ.d(0 as a function of the angle variable t for 
two different values of the scaled beam size Ox,s- When the beam size is larger filling more of 
the aperture, the intensity increases at all observation angles. While the absolute change in the 
maximum intensity is larger with larger beam size, the relative change in the minimum intensity 
(at ? = 0) is larger than that in the maximum intensity. 

This was exploited in the KEK experiments [1] where the beam size was determined from 
the ratio of the minimum to the maximum of the angular spectrum and we apply the same 
technique here. The top plot in Figure [10] shows the ratio of the minimum to maximum of the 
horizontal and vertical polarization intensities, d^S^ /dQ.d(0 and d^S^ /dQ.d(0 respectively, as a 
function of the scaled horizontal beam size o^^s- We mention an important point here: the ratio 
is independent of the energy, bunch intensity etc but depends only on the scaled parameters 
[g, u,t) and scaled beam parameters {Ox.s, <7v...s,-'co,i-,yo,5)- Hence as long as the scaled variables 
have the same values, these ratios will be the same for the Tevatron, LHC etc. The figure shows 
the data points calculated from the expressions in Equations (l55l) . (|6TI) and (l62l) as well as a 
quadratic fit through these points. The ratio for the horizontal polarization increases quadrati- 
cally with the horizontal beam size but the vertical polarization is insensitive to the change in 
horizontal beam size. Since the fit to the horizontal spectrum is quadratic, a 1% increase in the 




(62) 
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horizontal beam size results in a 2% increase in the minimum to maximum ratio. This suggests 
that if intensity differences at the level of 2% can be resolved, then beam size changes at the 
level of 1% can be detected. The bottom plot in Figure [10] shows the ratio for the horizontal 
polarization at three frequencies: u = 0.5, 1. The ratio increases quadratically with Ox at all 
frequencies. Furthermore, the ratio increases with u showing that the sensitivity of the beam 
size measurement increases with frequency. 

We consider now the sensitivity of the angular spectrum to changes in the beam offsets. To 
be useful for diagnostics, the spectrum should be sensitive to changes in offset which are frac- 
tions of a beam size. We find that when the offsets are in the range (0 — 1 ) x C7, the maximum of 
the angular spectrum does not change significantly. However the minimum of the spectrum at 
t = does change rapidly. So we examine instead the relative change in the minimum intensity 
as the offset changes. Figure [TT] shows the relative change of d^S^{t = 0)/dQ.d(O to changes in 
the beam offset shown in units of the beam size for three values of the frequency u. We observe 
that as the horizontal offset increases from zero to one times the beam size, the minimum hor- 
izontal polarization doubles in value. The relative change in the minimum is nearly the same 
for the three values of u, so at these frequencies the relative change in the minimum is not very 
sensitive to frequency. The bottom plot in Figure [TT] shows the ratio of the minimum to the 
maximum of the horizontal polarization. This ratio also increases quadratically with the offset 
and more importantly is sensitive to the frequency, increasing at higher frequencies. These re- 
sults show that if changes in the minimum intensity at the level of a few % can be measured, 
then changes in beam offsets of fractions of a beam size can be detected. 



5 Photon yield 

We start by calculating the differential spectrum for a single particle that is offset from the 
center of the hole. This will generalize the results in Section [3.1.11 We will then use this result 
to calculate the photon yield from a bunch and its dependence on beam and target parameters. 

The differential spectrum for a single particle is found by integrating the single particle 
differential angular spectrum over the solid angle. 



dW f d^W 2 fy f'^^ d^W t 

-l—\part = / ,„ , SmepdGpd^p=^ dt d^p———\part^=== 

da J dQ.d(Opart J Jo Jo dQ.d(o ^\-f-jyi 

Substituting from Equation ([48l) and integrating over the trigonometric functions we find 



\d<^pY,^^n = m'J'-f)f + {^plC-fW^^^ (64) 

where bm,n is the Kronecker delta and we have used U,n-m = { — ^)'^U,n,m- Note that the de- 
pendence on X has disappeared after this integration over ^p. After some simplifications we 
find 



2 



-^\p-rAS.U,bs)- — 



'2 , ^rnjJ_ 2j2^ 
ubj '"^ 



(65) 
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Figure 10: Top: The ratio of the minimum/maximum horizontal polarization intensity and verti- 
cal polarization intensity vs the scaled horizontal beam size Ox s at constant w = 1 , ^ = 1 .3, Oy,s = 
OA.xqs = 0.01 = jo,i- The horizontal polarization intensity increases quadratically with the 
beam size as shown by the quadratic fit while the vertical polarization intensity is insensitive 
to the horizontal beam size changes. Bottom: The minimum to maximum of the horizontal 
polarization intensity for two values of u : 0.5, 1. At m = 0.5, the ratio is smaller which shows 
that the minimum intensity is less sensitive to the beam size at the lower frequency. However 
the ratio still increases quadratically with the beam size. 
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Figure 11: Top: The change in the minimum of the horizontal polarization intensity in % as a 
function of the the horizontal beam offset, shown in fractions of the beam size for three values of 
u. This change is not sensitive to the frequencies in this range. Bottom: Ratio of the minimum to 
the maximum horizontal polarization vs the offset for the same values of u. This ratio increases 
with the frequencies. 
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where we have defined the integrated functions 

fy t 1 

Um,mig,u,r)= dt o ; I ,2^2 ^'"."'(g'^'0 (66) 

The differential spectrum from a bunch is found from 

dW f f dW N f _ dW 

{-T-)\bunch{g,u,r) = dbb dXp{b,x)-T- = „ / dhshsps{hs)— (67) 

a(a J J dco part 2;rc7x,.s<7y,j J dco pan 

where we have integrated over the scaled density 

Psibs) = Jdxpsibs.x) (68) 



From the expression for the density in Equation (1591) it follows that we can write 

Ps{bs,x) = exp[--(-^ + -^)]exp[-Z>,2c7+,,]exp[-Z>?a_,,cos2;f] 

xexp — y^cos;t: exp — 7^sm;t: (69) 
0"x,. fx,. 

Here we have defined the scaled beam parameters 

Expanding the exponentials of the trigonometric terms in Bessel functions, we have 
exp[-Z,?cT_,cos2z + ^cosz + ^sin;t] = £ £ I(-l)^'^P.A/;(^?a_,)4(^) 

^x,s ^x,s j>0k>0l>0 

bsJO s 

x//( — j^)cos27Xcos^;t:cos[/(;r/2-;i:)] 

where 

Dj = \ for j = 0, Dj = 2 forj>l 
After integrating over the angle x, it follows that 

p,ib,) = ^exp[-l(^ + ^)]G(Z..) (71) 



G{bs) = exp[-Z?^o-+,.,j 



4^ [bs ,s)h ( ^^2- ( —2- 



+ I £(-l)^+'Dp2//;(^^CT_,)/2K^) |^2|;-/|/2|,-Z|(^) 



A'o/[j=0=/l />0 >',^ 

bs^o 

+D2{j+l)h{j+l){—2 



(72) 
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If the offsets vanish, i.e. xo = = yo, this simplifies to 



lim ps = 2;rexp[-^^^+,,]/o(Z7^c7_.,) 
Define the following function obtained after integrating over bs 



(73) 



Fmiu; (7x,s: C7y,i, X^i, Jo,. 







/3 



(74) 



In the limit that the offsets vanish, this simplifies to 

■1 



lim F„ = 4/ JZ7.Z7,exp[-Z.2a+,,]/2(Z.2a_,,)[(4(^))2 + (^)2/2(^)] (75) 

The differential spectrum from a bunch is obtained from 

2 



(76) 



The photon yield into a bandwidth Act) is 



m>l 




(77) 



This is the photon yield from a single bunch per turn over the full 4n solid angle, the number 
of photons intercepted by the detector will be reduced by the acceptance of the detector. 

Figure [I2l shows the photon yield (plotted on a logarithmic scale) from a single bunch per 
turn into a 1% bandwidth as a function of the scaled frequency u = co/cOc for three values of 
the scaled beam size Ox.s- We have set the offsets xo.yo to zero for this calculation. The photon 
count is about 10^ photons per bunch per turn at w = 1 which should yield a detectable signal. 
The dependence of the photon yield on frequency can be fit to exponential curves with two 
different exponents below and above u = I. For example, for Ox,s = 0.125 we find that 



ANph ~ exp[-1.28M], u<l 
~ exp[— 1.82m], u>l 



(78) 



These exponents are not very sensitive to the beam size cr^ 5, e.g the exponents below and above 
M = 1 are (—1.29m, —1.85m) for a^,.? = 0.1. 

The photon yield does not depend sensitively on the beam size at low frequencies but at 
higher frequencies, the sensitivity to beam size increases. 

The scaled frequency u can be converted to a physical frequency by making a choice of the 



inner radius a/„ of the hole. If we assume o 



x,s 



0.125 



Oy s or ain 



So and o ^ 0.4mm at 



the CO location in the Tevatron, then the critical wavelength Xc corresponding to the critical 
frequency cOc is = 2;ra,„/7= 19/im. This is in the infra-red range. Detection in the optical 
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Figure 12: Photon yield (on a log scale) from a single bunch per turn into a 1% bandwidth as 
a function of the scaled frequency u for three values of the scaled beam size Ox^s- In all cases 
= <yx,s- Tevatron bunch intensity and energy were used in this calculation. The photon 
yield is not very sensitive to the beam size at low frequencies but the relative difference in yield 
increases at higher frequencies. 



le+08 



le+07 



o 
o 
-a 




Figure 13: Photon yield (on a log scale) from a single bunch per turn into a 1% bandwidth as a 
function of the scaled frequency u for two values of the target size g. Here the scaled beam size 



is chosen as Or 



0.125 = a 
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range seems to be ruled out since we find that at m = 19 or A = 1/im, the photon count is 
about 10^^/bunch/tum or practically zero. We can increase the photon yield by increasing 
the target size to say g = 2.5, the optimal value found earlier. Figure [13] shows the yields 
for ^ = 1.3,2.5. At low frequencies m < 1.5 the difference is significant but not so at higher 
frequencies. There is therefore no advantage to be gained with a larger target size if we choose 
to operate at frequencies above u> 2. 

We can use the photon yield to choose the frequency at which to observe the ODR. In 
practice the minimum photon yield will be determined by the efficiency and resolution of the 
camera and the level of the background synchrotron radiation which should be below the ODR 
photon yield. If for example we set the minimum photon yield to be 10^ photons/bunch/tum 
into a 1% bandwidth, then the highest scaled frequency from the plots above is m = 3.5. For the 
Tevatron, this implies that observable wavelength has to be at or above 5.4/im. If we integrate 
the signal over all 36 bunches and the bandwidth is greater than 1%, then the wavelength could 
be reduced some more. However it will be in the few microns range and not in the optical range. 



6 OTR vs ODR spectrum 

OTR is generated when the beam goes through the material of the target. The expressions for 
the OTR spectrum can be found from the ODR spectrum by taking the limit a„, — > 0. With 
OTR we cannot define a critical frequency (Oc and OTR observations show that there is no 
frequency at which the frequency peaks. This is a qualitative difference from the ODR spectrum. 
Similarly we cannot also define universal expressions for the OTR spectrum which depend only 
on dimensionless parameters. 

In this section we will briefly discuss the OTR single particle spectrum as a limiting case 
of the ODR spectrum. First consider the single particle at the center of the OTR target. From 
Equation (fTSi) it follows that the differential angular spectrum in the far-field is given by 

1 Q^y^ — — — 

= - (-1 ^^2-)2 \}^h{^ci)Kx{ad) + aJi{kd)KQ{ad)f (79) 

Here a is the radius of the target, k — k'&mQp, a = kly. 

Figure [14] shows the dependence of the single particle OTR and ODR spectrum on the 
angular variable t. For the parameters chosen here (Uout = l-3a/,,, u = 1), the maximum OTR 
intensity is about 8 times larger than the ODR intensity. Of this, a factor of 2.45 is due to the 
different areas of the material in the target for OTR and ODR. Furthermore, the OTR spectrum 
peaks at a smaller angle ?~1^0p~l/7 while the ODR spectrum peaks at ? ^ 1.6. This 
figure shows that ODR is beamed in a broader cone and at a larger angle than OTR. For the 
single particle spectrum with the particle offset from the center. Equation (15T] ) is still applicable 
for the differential angular spectrum except that U,„n defined in Equation (148]) simplifies to 

Um,n = a^[kJm-i(ka)Km{aa) + aJm{ka)Km-i{aa)] 

x[kJn-i(ka)K„{aa) + aJ„{ka)Kn-i{aa)] (80) 
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Figure 14: Comparison of ODR and OTR intensities for a single particle centered on the target. 
The outer radius aout is assumed to be the same in both cases. 

The OTR bunch spectrum can be similarly found by taking the limit a;„ — >^ in the appro- 
priate expressions above. 

7 Experimental Issues 

The devices and experimental conditions needed to observe ODR in the Tevatron requires a sep- 
arate detailed discussion. Here we will only mention some of the issues. Synchrotron radiation 
from the upstream dipoles hitting the target is an important source of background and needs to 
be mitigated. Preliminary calculations jSll show that the level of background at a target near 
the CO point in the Tevatron is less than the anticipated ODR flux. A mask placed upstream of 
the target may help to reduce this background to acceptable levels. The wavelength at which to 
observe the ODR depends on several competing factors. At longer wavelengths the ODR flux 
is higher but far infra-red detection is less sensitive, slower and complicated by other matters 
such as choice of windows which are sufficiently transparent at longer wavelengths. The syn- 
chrotron radiation background also increases at longer wavelengths. A satisfactory compromise 
might be in the vicinity of 5/im. Given the likely speed with which the ODR images will be 
acquired, it is unlikely that bunch by bunch and turn by turn imaging will be possible with the 
ODR monitor. Averaging over turns will most likely be necessary but it should be possible to 
update images on the time scale of seconds. Calibration of the ODR measurements requires 
measurements by other beam imaging devices nearby. At CO, the synchrotron light monitor is 
relatively close and would be suitable for calibration of the ODR monitor. If a circular target 
is used, it would likely be made of two semi-circular halves which will be moved in towards 
the beam on separate stepper motors. These two halves will need to be precisely aligned and 
their positions measured with respect to the beam. It is very likely that due to the nature of the 
helical orbits, the target in the Tevatron will only be suitable for one beam. We would choose 
the proton beam for ease of availability and greater intensity. It would be preferable but not 



27 



essential to choose the longitudinal location of the target so that both beams are not present 
simultaneously to avoid any effects from their parasitic interaction. Independently of this, the 
circular target may not be suitable if the helical orbits are separated by several beam sizes at the 
target. In this case the desired beam may have the required separation from the target but the 
opposing beam will not be far enough from the target. This could be avoided by having the two 
halves of the target separated by a gap. Instead of semi-circular foils, rectangular foils on either 
side may be preferable in this case. 

We envisage that if initial measurements are successful, the ODR monitor could be used a 
passive device monitoring beam parameters and their changes during the length of a store. Here 
other operational challenges will arise. For example, beam motion which changes the position 
relative to the target will need to be included in the automated ODR measurement. Wake fields 
due to the target and heating of the target by beam induced fields are likely to be negligible but 
need to be considered. Some of these same issues have arisen and been resolved with the use 
of the pick-off mirrors for the synchrotron light mirror and will benefit from that experience. A 
detailed account of these and other relevant issues will appear elsewhere. 

8 Conclusions 

Our main concern here was the far field ODR spectrum from a round hole in a collider and 
specifically the Tevatron. We list the main conclusions 

• Existence of a critical frequency and universal curves. 

There is a critical frequency cOc associated with ODR at which the angular spectrum 
intensity peaks. There is no such frequency with OTR This is a well known phenomena. 
However we have also shown, something not previously recognized, that the spectrum 
for a round hole depends only on dimensionless parameters. Hence the results seen in 
Figures |7] to [HI especially the sensitivities, are universally applicable to all machines 
when these dimensionless parameters have the same values. 

• Sensitivity to beam sizes. 

The ratio of the minimum to maximum intensity is very sensitive to the beam size as 
seen in Figure \W\ The sensitivity increases with observation frequency. The horizontal 
polarization is sensitive to the horizontal beam size and similarly for the vertical plane. 
This differs from the dependence with rectangular slits and straight edges. 

• Sensitivity to beam offsets 

We found that the minimum of the angular spectrum along the direction of specular re- 
flection is most sensitive to changes in the beam offset. The ratio of the minimum to the 
maximum can also be used, this has the added advantage of higher sensitivity at increas- 
ing frequencies. Again, we found very good sensitivity (seen in Figure [TT]) to changes in 
beam offset of fractions of a beam size. 
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• Photon yield. 

We calculated the photon yield from a single Tevatron bunch at 980 GeV. The photon 
yield decreases exponentially fast with frequency. Assuming an 8a separation between 
the beam and the target, the critical wavelength is 19/im. The calculation predicts that 
photon yields of ~ 10^ photons/bunch/tum into a 1% bandwidth will be obtained at about 
a 5/im wavelength. Detection at optical frequencies does not seem feasible. It is clear 
that the lower wavelength limit for observable ODR signals in the far-field is at a few 
microns, in the infra-red regime. 

In this paper we did not consider the near-field spectrum in much detail except briefly in Sec- 
tion 13.1.21 That discussion as applicable to the spectrum from a bunch and a detailed discussion 
of experimental details will appear separately 
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